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METHOD FOR SELECTING REPRESENTATIVE ENDMEMBER 
COMPONENTS FROM SPECTRAL DATA 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The subject matter of the present invention generally involves spectral image 
analysis. More particularly, the present invention involves processes for determining 
representative scene components, or endmembers, from spectral imagery. 

Description of the Related Art 

It has long been recognized that at spatial scales typical of remotely sensed data, 
the surface corresponding to a given pixel is rarely composed of a single physical 
component, but rather a combination of "sub-pixel" constituents. This mixed-pixel 
concept was developed into a quantitative model for interpretation of spectral images, 
Spectral Mixture Analysis ("SMA"). SMA involves an interactive selection and 
refinement of key components, or endmembers, to achieve an acceptable solution. The 
SMA model is based on the assumption that remotely sensed spectral measurements are 
mixed signatures, which vary across the scene as the relative proportion of each 
component changes. SMA can therefore be used to translate spectral imagery into 
quantitative maps also referred to as abundance maps, showing the percent cover of each 
spectrally distinct scene component. The number of endmembers resolvable with SMA 
is dependent on the complexity of the actual surface, the spatial resolution of the imagery, 
and the spectral resolution of the data. However, even in hyperspectral imagery of 
complex regions, SMA is most robust when a few (<10) spectrally distinct components 
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are modeled. For instance, a forested scene would be modeled as a combination of green 
vegetation (e.g. leaves, grass), non-photosynthetic vegetation (NPV, e.g. bark, wood, dry 
leaves), soils, and shade. This last component, shade, is unique to SMA. The shade 
endmember accounts for spectral variations in illumination caused by either topography 
5 or sub-pixel surface texture. With SMA, a tree or group of trees can therefore be modeled 
not by a single "tree" signature, but instead by a combination of green vegetation, NPV, 
and shade. The selection of these types of endmembers generally requires at least one 
expert user's input to be successful. 

Alternatively, complex statistical and mathematical approaches have been 

10 employed in order to perform endmember selection, but none have been shown to 
provide routine and reliable results. These approaches utilize purely mathematical 
transformations of the spectral data, without regard for what physical information the 
data actually represents, e.g., vegetation, NPV, soil, shade, buildings, roads, etc. The 
endmembers determined through these processes are physically meaningless and not 

1 5 comparable to those derived by expert users. As such, their utility is limited and at best 
require subsequent subjective interpretation by an expert. 

Current endmember selection methods are thus either manual with expert users 
doing the processing or automatic using signal processing, but with less reliable results 
that still require expert user interpretation. 

20 
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SUMMARY OF THE INVENTION 

Summary of the Problem 

Currently, the process of determining the number and nature of these spectrally 
distinct components (or "endmembers") using either the manual or statistical selection 
process is both subjective and labor intensive. Because of this, SMA-based 
classification, although clearly of significant use in many applications, has largely been 
limited to labor intensive research projects. In such studies, endmembers are typically 
selected manually by experienced analysts. ^Results from this process are subjective and 
vary from one investigator to another. Further, as the complexity of data sets increase, so 
does the time required to perform SMA as does the likelihood of greater error. 
Consequently, there is a need in the art for a faster, more objective, more cost-effective, 
user-friendly process for selecting representative spectral endmembers for a scene. As 
more and more data become available from, e.g., commercial airborne systems and earth 
observational satellites which include hyperspectral sensors, there is an increasing need 
for a process that aids both experienced and novice end-users to work more efficiently, 
significantly reduces costs, and enhances the utility of spectral imagery. 

Summary of the Solution 

The invention described herein is an objective and automated process for 
determining spectral endmembers and transforming Spectral Mixture Analysis (SMA) 
from a widely used research technique into a user-friendly tool that can support the needs 
of all types of remote sensing. The process extracts endmembers from a spectral dataset 
using a knowledge-based approach. The process identifies a series of starting spectra that 
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are consistent with a scene and its environment. The process then finds endmembers 
iteratively, selecting each new endmember based on a combination of physically and 
statistically-based tests. The tests combine spectral and spatial criteria and decision trees 
to ensure that the resulting endmembers are representative of the scene. 
5 The invention described herein achieves a significant advancement in the 

application of SMA. As discussed herein, the process of the present invention 
successfully implements a number of the manual decision making steps that an analyst 
performs in doing mixture analysis. However, these steps are ordered in a way that leads 
to far fewer false steps, fewer refinements of endmember selection, and, minimizes the 

10 need to remove entire endmember sets from an analysis sequence. Further, the process of 
the present invention provides users with sophisticated insight into mixture model 
methods. This is particularly important when working with, for example, complex scenes 
or on areas where specific objectives must be met. 

As a result of the processes described herein, the time required to exploit spectral 

1 5 image data is dramatically reduced. The process makes information from spectral 
imagery readily available, and expedites critical decision making. The process 
automatically determines the subpixel proportions of materials in spectral datasets and 
transforms these data into readily interpretable information layers that can be imported 
into, e.g., a geographical information system (GIS) database. 

20 Applications of spectral imaging data can support a diverse range of scientific and 

business operations. Among these, spectral data offers enormous potential for routinely 
creating geologic and soil maps, environmental monitoring, aiding forestry and 
agriculture production, land use/land cover mapping, and various applications in support 
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of civil, military, law enforcement, and intelligence communities. For example, high 
resolution spectral imaging is used to help winegrowers harvest their crops with more 
precision. By identifying plant health and maturity, vintners can establish optimum 
harvest schedules that ultimately yield high quality wine. The same methodology holds 
5 true for using spectral imaging data to non-destructively detect bruised fruit in large 
orchards, and also aiding law enforcement officials in detecting marijuana plants from 
airborne platforms. 

More particularly, an embodiment of the present invention provides a process for 
determining at least one candidate spectral endmember that represents a group of N 

10 spectra. The process of this embodiment comprises: building an initial endmember set 
composed of at least a first and a second spectrum representing at least a first and a 
second spectral characteristic expected to be in the group of N spectra; unmixing the N 
spectra within the group of N spectra to determine what portion, if any, of each of the N 
spectra match at least one of the at least a first and second spectrum; calculating a first 

1 5 metric value for each of the N spectra, wherein the first metric value accounts for a 
remaining portion of each of the N spectra not matching the at least a first and second 
spectrum; defining a metric value range, wherein the N spectra having first metric values 
within the metric value range are defined as M spectra; ordering the M spectra from 
highest first metric value to lowest first metric value; comparing each of the M spectra, 

20 beginning with the M spectra having the highest first metric value, to each of the N 

spectra, to determine the frequency with which each of the M spectra occurs within the N 
spectra; and calculating a second metric value for each of the M spectra, wherein the 
second metric value combines the frequency of occurrence of each of the M spectra 
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within the N spectra with a first metric value for each of the M spectra, wherein the M 
spectra having the largest second metric value is the at least one candidate endmember. 

A further embodiment of the present invention provides a process for determining 
at least one candidate spectral endmember that represents a group of N spectra. The 
5 process of this embodiment comprises: building an initial endmember set composed of at 
least a first and a second spectrum representing at least a first and a second spectral 
characteristic expected to be in the group of N spectra; unmixing the N spectra within the 
group of spectra to determine what portion, if any, of each of the N spectra match at least 
one of the at least a first and a second spectrum; defining an error value for each of the N 

10 spectra, wherein the error value is the portion of each of the N spectra that does not match 
a combination of the at least a first and a second spectrum; comparing the error value for 
each of the N spectra to a predetermined error value range, wherein spectra having error 
values within the predetermined error value range are defined as M spectra; ordering the 
M spectra from highest error value to lowest error value; comparing each of the M 

1 5 spectra, beginning with the M spectra having the highest error value, to each of the N 

spectra, to determine the frequency with which each of the M spectra occurs within the N 
spectra; and calculating a metric for each of the M spectra, wherein the metric combines 
the frequency of occurrence of each of the M spectra within the N spectra with an error 
value for each of the M spectra, wherein the M spectra having the largest metric is the at 

20 least one candidate endmember. 

A still further embodiment of the present invention provides a process for 
determining at least one candidate spectral endmember within a scene having N pixels 
using scene spectral data. This embodiment of the present invention comprises: building 
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a first endmember, wherein the first endmember represents a first spectral characteristic 
expected to be in the scene; building a second endmember for the scene, wherein the 
second endmember represents a second spectral characteristic expected to be in the scene; 
unmixing the N pixels in the scene to determine what portions of each of the N pixels 
5 match at least one of the first and second endmembers; defining a remaining portion for 
each of the N pixels not matching a combination of the first and second endmembers as 
an error value, wherein each error value corresponds to a residual spectra and the residual 
spectra in combination form a residual spectrum of the N pixels; calculating a root mean 
square (RMS) error for the N pixels by combining the error values for the residual 

1 0 spectra; determining an acceptable range of deviation from the mean RMS error; 

comparing each of the RMS error values for each of the N pixels to the acceptable range 
of deviation from the mean RMS error and keeping the M pixels that are within the 
acceptable range of deviation; ordering the M pixels from highest RMS error value to 
lowest RMS error value; comparing the corresponding residual spectra of the M pixels, 

1 5 beginning with the M pixel having the highest RMS error value, to the residual spectrum 
comprising the residual spectrum for the N pixels, to determine the frequency with which 
each of the corresponding residual spectra of the M pixels occurs within the residual 
spectra for the N pixels; and calculating a weighting factor for each of the M pixels, 
wherein the M pixel having the largest weighting factor is the at least one candidate 

20 endmember. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings: 

Figure 1 is a flowchart for finding scene reflectance spectra endmembers according to an 
embodiment of the present invention; 

Figure 2 is a flowchart for initiating an endmember sequence according to an 
embodiment of the present invention; 

Figure 3 is a flowchart for initiating a display sequence according to an embodiment of 
the present invention; 

Figure 4 is a flowchart for initiating an unmixing sequence according to an embodiment 
of the present invention; 

Figure 5 is a flowchart for initiating a sequence for finding endmembers according to an 
embodiment of the present invention; and 

Figure 6 is a flowchart for initiating a sequence for displaying status histograms 
according to an embodiment of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS OF THE PRESENT 

INVENTION 

The following embodiments, though representing the best known mode of 
practicing the present invention, are not intended to be limited thereto. Based on the 
knowledge of one skilled in the art, variations of the following embodiments are intended 
to be part of the scope of the present invention. 

The flowchart of Figures 1-6 sets forth the general process steps for determining 
the land cover types, e.g, green vegetation (leaves, grass, etc.), non-photosynthetic 
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vegetation ("NPV") (e^, bark, wood, dry leaves), soils, snow, ice, roads, buildings, and 
shade, within a particular scene according to an embodiment of the present invention. 
Image cube reflectance, emission, or radiance data (hereafter "reflectance data") 
indicative of a scene is read into a processor S10. The scene boundaries are defined by, 
for example, a set of ground-based longitudinal and latitudinal coordinates. The 
reflectance data results from imaging performed by any of a number of sources, such as, 
satellites, aircraft, weather balloons, in situ or ground based sensors, or the like. One 
skilled in the art recognizes that the current invention may be used with all types of 
spectral data, ( e.g., multispectral, hyperspectral, ultraspectral) with results varying in 
level of detail depending on spectral and spatial resolution and spectral range. 

As described above, no matter the specifics of the datasets, a spectral 
measurement in a given pixel is at some scale mixed. For example, even with state of the 
art remote imagery, in order to identify and classify the land cover types that are 
represented in the selected scene, it is necessary to determine spectrally distinct 
components from the combination components received by the imaging system. These 
spectrally distinct components for a particular scene are referred to as endmembers. 

Again, referring to Figure 1, in an embodiment of the present invention, after the 
reflectance data for a particular scene is read into the processor S10, pre-existing 
endmember data for that particular scene is also read into the processor, if available S30. 
The process of the current invention stores any previously generated endmember data 
according to location within the scene and the spectrum for each endmember as a 
function of wavelength. Consequently, after the initial reflectance data is read into the 
processor S10, the user is queried as to whether the user chooses to load in an existing 
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endmember file for the selected scene S20. If the user chooses to load in an existing 
endmember file, the file is read into the processor, in for example, ASCII format S30. 
The endmember data is subjected to a programming subroutine in order to format the data 
into a recognizable and easily manipulatable endmember data structure S200. 
5 Alternatively, if there is no existing endmember file for the selected scene or the 

user chooses not to load the existing endmember file, the user initiates an endmember 
sequence S100 in order to create a new endmember structure. Referring to Figure 2, 
after initiating the endmembers sequence SI 00, a shade spectrum, ejj., 1% percent 
reflectance spectrum, i.e., unchanging with wavelength, is constructed at the channels for 

1 0 the scene reflectance data SI 20. In the embodiments described herein, a user assumes 
that a scene's reflectance data is representative of at least two components, shade or 
shadow (hereafter "shade") and vegetation. A shade spectrum is constructed in order to 
take into account and adjust for the presence of inhibiting shade components and can be 
considered the first endmember. The second endmember is determined from a reference 

1 5 vegetation spectrum that has been read into the processor and resampled using linear 
interpolation to match the wavelengths of the reflectance scene data set S140. The 
reference vegetation spectrum is selected from a library of vegetation spectra in an 
attempt to find a match between known vegetation spectra and scene reflectance data. 
These preliminary endmembers, i.e., shade and vegetation, are included in the 

20 preliminary endmember structure S160 that is outputted at S180 and returning to Figure 
1, is used in a display sequence S300. Alternative embodiments, use various other initial 
endmember sets, consisting of at least two of the materials expected to be in the scene, 
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based on known data and scene properties. Other materials might include soils, NPV 
(such as bark, wood and dry leaves), snow, etc., as listed previously. 

To this point, the process steps described above and the results thereof have been 
invisible to the user. In order for the user to view results and change parameters 
5 accordingly, an embodiment of the present invention facilitates a display of the results. 
Steps for initiating a display sequence according to an embodiment of the present 
invention are shown in Figure 3. The preliminary endmember structure and the scene 
reflectance data are inputted into a display sequence ( e.g. , computer program) S3 10. In 
order to provide visual aids to the user, an embodiment of the invention provides both a 

10 true color image from specified bands in the scene reflectance data S320 and a 

normalized difference vegetation index ("NDVI") from specified bands in the scene 
reflectance data S330. NDVI is the traditional vegetation index used by researchers for 
extracting vegetation abundance from remotely sensed data. It divides the difference 
between reflectance values in the visible red and near-infrared wavelengths by the overall 

15 reflectance in those wavelengths to give an estimate of green vegetation abundance. In 
essence, the algorithm isolates the dramatic increase in reflectance over the visible red to 
near infrared wavelengths, and normalizes it by dividing by the overall brightness of each 
pixel in those wavelengths. Specifically NDVI is: 

NDVI = IR-red 

20 IR+red 

The result of this algorithm is a single band data set, ranging from -1 to 1, with values 
corresponding to photosynthetic vegetation abundance. True color and NDVI images are 
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merely exemplary of the types of images which may be presented to a user during the 
display sequence. One skilled in the art recognizes the various images which fall within 
the scope of the invention. 

Next, similar to steps recited with reference to Figure 1, the display sequence 
5 queries whether or not there is an existing endmember file, or if the sequence must 

assume that there is no pre-existing information and start from scratch to determine initial 
endmembers for the scene S340. In the case where there is no pre-existing endmember 
information, in order to identify initial endmembers, spectral matching algorithms such as 
Spectral Angle Mapper ("SAM"), or the like, are used to find the closest image spectrum 

10 in the scene reflectance data to reference vegetation spectra S350. In order to ascertain 
the closest spectral signature to specified bands of the scene reflectance data, SAM 
measures the similarity between an unknown spectrum and the reference spectrum. SAM 
treats two spectra as two vectors in n-dimension space, where the dimensionality equals 
the number of bands. The similarity of two spectra ( e.g. , vectors) will be based on the 

15 angle of separation between the two vectors. The smaller angles represent closer matches 
to the reference spectrum. After running the scene reflectance data through a matching 
program, such as SAM, the outputted endmember structure from SI 80 is updated to 
include the new vegetation spectrum data S360 and a spectral unmixing sequence is 
initiated S370. 

20 Referring to Figure 4, after initiating the unmixing sequence S370, endmembers 

are compiled from the updated endmember structure S372 and spectral mixture analysis 
is performed on the compiled endmembers S374. The technique of pixel unmixing as it 
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is generally applied assumes a linear model to separate the primary scene components, 
termed endmembers. Alternatively, non-linear models are applied. The technique relies 
on the analyst to identify the spectra of the endmembers. These spectra are then typically 
used in a least-squares inversion to estimate the contribution of each endmember for each 
5 pixel by simultaneously solving the following equations for each wavelength b: 

m 

1=1 

. m 

S /, = 1 o 

i = 1 

Rb = input pixel spectrum at wavelength b 
fi = contribution or fractional abundance of the ith endmember 
1 0 rib = spectrum of the ith endmember 
Eb = a residual error at wavelength b 

The model is that Rb is equal to the sum over the number of endmember of plus 
1 5 Et. Variable is the endmember; f is the abundance of the end member i\ Rb is the pixel 
spectrum for each end member at band b\ and Eb is the residual error at each band b. A 
second constraint, that the sum of fi will equal unity (1.0), may also be applied. The 
results of applying this analysis are a series of images, one for each endmember, 
containing the percentage, or fractional abundance, of the endmember in each pixel. The 
20 collection of abundance images for all endmembers are often referred to as fraction 
cubes. 

Depending on the results of the SMA step, a shade mask is built and applied to 
the unmixed pixels, ej^, all pixels with greater than a user selection percentage of shade, 
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e.g. , 80% shade are masked out S376. Next a root mean square ("RMS") error image is 
generated and it replaces the NDVI image that is viewable by the user S378. The RMS 
error image is an image displaying the square root of the sum of the squares of the error 
as a function of wavelength for each pixel. In an embodiment of the present invention, 
5 the RMS error image is an intensity image where bright pixels are areas that are not well 
defined by the current choice of endmembers, i.e., shade and initial vegetation. The 
endmember spectra are also plotted in a graph, wherein, for example, the x-axis is 
wavelength and the y-axis is reflectance S380. At any point after initiating the display 
sequence S300, the created and updated images and plots are viewable through display 
10 S400. 

Using the display, the user is able to visualize the images and plots in order to 
determine how to vary the parameters, e^, shade threshold, SAM threshold, RMS 
minimum and maximum (described below), in order to find the next endmember 
spectrum. The shade threshold parameter adjusts how many dark pixels to allow in the 

1 5 analysis. The value is based on the abundance value of the "shade" endmember. Areas in 
the "shade" abundance image that are greater than or equal to this value are masked out 
and are not considered as possible endmembers. The display S400 allows the user to 
view how the shade threshold constraint affects the spatial analysis of the image by 
displaying a Shade Mask window. Changing the SAM threshold parameter affects how 

20 many pixels in the image are flagged as spectrally similar to the spectrum of the new 
candidate endmember. The display S400 allows the user to view how this constraint 
affects the spatial analysis of the image when the colored mask is displayed over the 
RMS error image. When starting a new image, this parameter is initially a small value 
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(e.g., the default starting value is 0.1). As each new consecutive endmember is found, this 
value is increased to loosen the matching constraint. Before accepting any new 
endmember, the user defines the fine balance between spectral uniqueness of each 
endmember versus the spatial extent of that endmember. 
5 The RMS minimum and maximum establish how many candidate pixels in the 

image are considered for the next endmember. These two values are the number of 
standard deviations from the mean RMS error value of the image. A user can find out 
how many pixels are in the established range and can examine the data to determine the 
best "RMS Min/Max Sigma" values (discussed further below), through the interactive 

1 0 display S400 via a "Plot RMS Error" button. A window with a histogram plot of the 

RMS errors of every pixel in the image is shown on the display. The histogram is labeled 
with a legend to indicate mean RMS error value of the image and indicate the number of 
candidate spectra listed. The candidate number is the sum of the RMS error values 
between the two appropriately marked min/max sigma values from the mean. 

15 After initiating the display and unmixing sequences as discussed above, the user 

initiates the find endmember sequence S500. Referring to Figure 5, the updated 
endmember structure is compared to any existing fraction cube data to determine if a 
corresponding unmixing has previously been performed, e^g. at S370-380, S505. If not, 
i.e., SMAflag=0, then the unmixing sequence is performed S515. Once fraction cube 

20 data that corresponds to the endmember structure is in hand, the shade mask is retrieved 
S520 and pixels in the RMS image that are flagged in the shade mask are masked out 
from the endmember sequence S525. The pixels remaining after being compared to the 
shade mask at S525 are subjected to an RMS calculation to determine the mean RMS 
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value at S530. The individual RMS calculations for each pixel are next compared to a 
user defined RMS range of minimum and maximum from the mean RMS values S530. 
The pixels having RMS values (RMSval) falling within this range are sorted in 
decreasing order, e.g., from highest to lowest error, and are stored according to the 
5 coordinates of the corresponding pixel S535. 

Next, the corresponding residual spectrum (K x ) for each RMSval is retrieved 
S540. More particularly, the residual spectrum is the residual error as a function of 
wavelength for each pixel in the scene. Each K x , beginning with Kj, corresponding to the 
highest error, is compared to the residual spectrum for every pixel in the image to 

10 determine the number of pixels within the image that are similar to K\. Depending on 
the number of pixels in the image as well as processing limitations, the user selects a 
predetermined number of K x 's which are evaluated before the endmember determination 
sequence continues. For example, in S545, x=500 before a histogram sequence is 
initiated S560. Consequently, 500 RMSvals are evaluated before the histogram is 

1 5 constructed. The evaluation of each RMSval requires comparing each K x for each 

RMSval using a matching algorithm, e.g., SAM (as described above) with the residual 
spectrum for each pixel in the scene, i.e., SAMresidual(K x ) = SAM(RMSval(K x )) S550 
and counting the number of pixels in the scene containing the RMSval according to a pre- 
established threshold for the SAM value, i.e., SAMresidual(K x ) ^SAM threshold S555. 

20 Referring to Figure 6, when x=500, a new display window is built which includes 

a plot of a histogram of RMS values versus SAM pixel count, i.e., RMSval vs. 
RMSfrequency S565. Weighting factors are calculated as follows to determine the pixel 
coordinate for the most likely endmember candidate S570. The factors are defined 
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wherein RMSval(K x ) is the RMS error value for the corresponding pixel, RMSval(min) is 
the minimum RMS error value for the K x spectra, RMSval(max) is the maximum RMS 
error value for the K x spectra, RMSfreq(K x ) is the value corresponding to the number of 
scene pixels that match the residual spectrum for the corresponding RMS error value, 
RMSfreq(min) is the number of occurrences of the least frequently occurring residual 
spectrum from the scene residual spectrum, and RMSfreq(max) is the number of 
occurrences of the most frequently occurring residual spectrum from the scene residual 
spectrum. 

Normval(K x ) = RMSval(K x ) -RMSval(min) 



RMSval(max) - RMSval(min) 
Normfreq(K x ) = RMSfreq(K x ) -RMSfreq(min) 



RMSfreq(max) - RMSfreq(min) 
wgtfactor(K x ) = sqrt( Normval(K x ) * Normfreq(K x )) 

candidate = MAX( wgtfactor ) 



The above-cited equations, or the like, act as a metric to determine the candidate 
endmember with the highest error and the largest frequency within the scene. After 
plotting the histogram and determining the best candidate endmember based on the 
25 weight factors, a color mask is created on an image of the scene in order to highlight for 
the user the pixels containing the candidate endmember S575. Using this color mapped 
image, the user can determine based on, among other things, spatial coherence, whether 
this candidate endmember is acceptable S580. If the candidate endmember is not 
acceptable, the process ends S585 At this point, the user can for example, select the 
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candidate endmember with the next highest weight factor or adjust the SAM threshold, in 
order to find new candidate endmembers. If the candidate endmember is acceptable, it is 
added to the endmember structure already containing shade and initial vegetation 
endmembers S590. Finally, using this new endmember with the two previous 
5 endmembers, the scene reflectance data is subjected to the unmixing sequence S370-S380 
and the process ends S600. This process is iterated to find new endmembers. 

In an alternative embodiment of the present invention, the first new acceptable 
candidate endmember replaces one of the original two endmembers and the unmixing 
sequence is run on the scene reflectance data with the two remaining endmembers, to find 
1 0 a new third endmember. 

In an alternate embodiment of the present invention, the whole of the processes 
described herein are automated, wherein all user interactions, definitions, and decisions 
(including how many endmembers there are) are completely replaced with decision tree 
algorithms. In this automated embodiment, preferable ranges are iteratively selected by a 
1 5 processor based on pre-programmed decision trees. 

The embodiments described herein are intended to be exemplary, and while 
including and describing the best mode of practicing, are not intended to limit the 
invention. Those skilled in the art appreciate the multiple variations to the embodiments 
described herein which fall within the scope of the invention. 
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